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We study boson stars in Brans Dicke gravity and use them to illustrate some of the properties of 
three different mass definitions; the Schwarzschild mass, the Keplerian mass and the Tensor mass. 
We analyse the weak field limit of the solutions and show that only the Tensor mass leads to a 
physically reasonable definition of the binding energy. We examine numerically strong field uj = —1 
solutions and show how, in this extreme case, the three mass values and the conserved particle 
number behave as a function of the central boson field amplitude. The numerical studies imply 
- - - that for u> = —1, solutions with extremal Tensor mass also have extremal particle number. This is 

a property that a physically reasonable definition of the mass of a boson star must have, and we 
prove analytically that this is true for all values of cu. The analysis supports the conjecture that the 
Tensor mass uniquely describes the total energy of an asymptotically fiat solution in BD gravity. 
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' Boson stars (or Klein Gordon geons) are asymptotically flat self gravitating conflgurations of zero temperature scalar 
particles (bosons) . The boson field is described by a complex wave function 5* and the matter Lagrangian possesses a 
U{1) internal symmetry. This symmetry leads to the existence of a conserved charge N which is interpreted as the total 
number of bosons. The first boson star solutions were found by Kaup Q and independently by Ruffini & Bonazzola 
1^. These authors examined spherically symmetric boson stars in General Relativity (GR) for which the only self 
I interaction term in the matter Lagrangian is the boson mass. They found that the solutions exhibit qualitatively 
I , similar behaviour to those for neutron stars and white dwarfs: they form a single parameter family whose ADM 
^ mass Madm and particle number N vary smoothly with the parameter and have coinciding maxima and minima. 



However, the ADM masses of the stars were found to be extremely small, of order M^j/a, where Mpi := y/h/G i 
the Planck mass and a is the mass of the bosons. Colpi, Shapiro and Wasserman |^ studied GR boson stars whose 
matter Lagrangian includes a quartic self interaction term. They found that this term dominates the pressure of the 
stars and increases the ADM mass to order M^^/a^. 
5^ \ The stability of boson stars in GR has been studied by, amongst others, Kusmartsev, Mielke & Schunck Q- 

who used an approach based on catastrophe theory. To use this method one must identify two conserved charges that 
vary smoothly with some parameter which labels the solutions and whose extrema occur at the same values of this 
parameter. Then a plot of one charge against the other will show cusps at the extremal points and the first of these 
cusps marks the onset of dynamical instability. For GR boson stars, the appropriate conserved charges are N and 
Madm^ and Jetzer |^ has proved that these two quantities do have coinciding extrema both for pure boson stars and 
for mixed boson-fermion stars. 

There has been renewed interest recently in generalising GR by adding one or more scalar (dilaton) gravitational 
fields. These extra fields appear in the low energy limits of super string and super-gravity theories. The siniplest low 
energy string theory effective action is approximated by the uj = —1 action of Brans Dicke (BD) theory 0], which 
includes one dilaton field. This choice of BD parameter has been ruled out by solar system experiments that imply 
that oj > 500. However, these experiments do not place such strong limits on more general Scalar Tensor (ST) theories 
in which the coupling parameter to is allowed to vary. Thus it is conceivable that the gravitational interaction can 
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be described by a ST theory in which « — 1 at some early time when string theory is still valid, and that this 
parameter has increased to it present value as the Universe has evolved to its present state. 

In the original (Jordan frame) formulation of BD theory, Lee [|j has shown that Madm is not a good description 
of the total energy of a spacetime. In particular, it is not conserved by an isolated source that emits gravitational 
radiation. To obtain a conserved mass on must modify Madm by adding another term that involves the derivative 
of the dilaton field. 

Several authors have studied boson stars within the framework of ST gravity. Gunderson & Jensen j^] considered 
boson stars in BD theory, focusing on w = 6 solutions, and they showed that these solutions are qualitatively similar 
to those in GR. Tao & Xue ||l^ discuss solutions to a scale invariant variation of BD theory in which the mass term in 
the matter Lagrangian couples to the curvature via the dilaton. This theory necessarily violates the Weak Equivalence 
Principle and the authors found a conserved charge associated with both the boson and dilaton fields. Torres pd| ] 
examined boson stars in more general ST theories while Comer & Shinkai Jl^ used catastrophe theory to examine 
the stabihty of ST boson stars. 

When one considers a ST boson star in a cosmological setting, the evolution of the background (cosmological) 
value of the dilaton field may differ from the evolution of the dilaton field near to the centre of the boson star. This 
phenomena was originally described by Barrow as "gravitational memory" jist , although he considered the effect of 
cosmological evolution on the horizon area of a black hole. However, as pointed out by Torres, Liddle & Schunck [pT| , 
ST boson stars also offer a means of testing theoretically the degree to which the gravitational memory hypothesis 
is true, since boson stars are also highly relativistic objects and have the added advantage of being singularity free. 
These authors have also studied the gravitational evolution of a BD boson star in a cosmological background p^ ] 
and they found that a star that is stable at the current epoch remains stable as one traces its evolution backwards in 
cosmological time. 

In this paper we examine how one defines mass in ST gravity and use boson star solutions to illustrate the properties 
of three mass definitions: the Schwarzschild mass (which gives the ADM mass in the asymptotic limit), the Keplerian 
mass (that quantity that determines the geodesic motion of non-self gravitating test particles) and the Tensor mass 
(the conserved quantity identified by Lee ) . For the sake of simplicity we consider only BD theory. In the spirit of 
the original formulation of the theory we assume that the Jordan frame is the physical frame. We focus on = — 1 
strong field solutions since in this case the field equations are sufficiently different to the GR equations for strong field 
scalar-tensor effects to appear. Also, as mentioned above, this value of the coupling parameter may be relevant to the 
study of boson stars in the early Universe. However, our analytical results are valid for any value of to. 

The plan of the paper is as follows. In Section || we give the field equations and discuss the boundary conditions 



that must be imposed to obtain boson star solutions. In Section [II we discuss three definitions of quasi-local mass 
adapted to the symmetries of the solutions and give an expression for a generalisation of the ADM mass that is valid 
in BD theory. In Section |^ we examine the asymptotic form of the solutions and show that in the asymptotic limit, 
the three quasi-local masses each tend towards one of the generalised ADM masses. In Section ^ we describe how 
the physical characteristics of the solutions change under a rescaling of the dilaton field. In Section ^ we discus the 
weak field limit of the solutions and derive an expression for the fractional binding energy in this limit. In Section 



VII we describe some of the properties of the strong field solutions, focusing on th e to = —1 coupling and highlighting 



the differences between the three mass values we are considering. In Section VIII we prove that for all to, solutions of 
extremal particle number are those of extremal Tensor mass. Finally, in Section IX we give some concluding remarks. 
Throughout this paper we use the sign and index conventions of |l6| . 

II. STRUCTURE EQUATIONS AND BOUNDARY CONDITIONS 

We take as our starting point the action for Brans-Dicke gravity minimally coupled to a complex boson field 5*: 



S = 



1 

(7^ - uoV"(jNa(f) - 7 (V'"1'*V^* + Vl-V^**) - :;72 



(1) 



where a is the mass of the boson particles and we have chosen units in which c = 1. The dilaton field (j) is dimensionless 
and the product e'^G measures the strength of the gravitational coupling. Varying the metric leads to the field 
equations 

G^, = (1 + c^)V^0V,0 - (1 + c^/2)5^,V"0V,0 - V^V,.^ + 5^,00 + SttGc^T^, (2) 

where 
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(3) 



is the energy momentum tensor of the boson field. Varying the dilaton and boson fields gives the dilaton wave equation 

2wn(/) = n + ujV(l)V^(f> (4) 

and the boson wave equations 

2 2 

The action possesses a global U{1) symmetry which implies the existence of a conserved current 

Jf' ^ ^gf" -■^V^^^J*) , V^J^ = 0. (6) 

Note that has no explicit dependence on (j): eqn (||) is identical to the GR definition since the boson field is 
minimally coupled to TZ. For an asymptotically flat spacetime and time-like current, eqn gives a conserved 
(time-independent) charge 

7V = y" d^xy/hn^J^" (7) 

where the integral is taken over an arbitrary space-like hypersurface with time-like unit normal and induced metric 
hfiv = g^f + n^TLv In contrast with the solutions discussed in ||l^, there in no conserved charge associated with 
the dilaton field. Although we have written the dilaton terms on the right of eqn (||), so that they contribute to the 
total energy momentum tensor, we interpret 4> as an extra component of the gravitational field, not as an additional 
particle field. 

We consider static, spherically symmetric equilibrium solutions. Using the standard orthogonal {i, r, 0, ip} coordi- 
nate basis we write the line element as 



di' -f e^^dH' + K'de' -h K' sin" e dip'' 



(8) 



where ly and A are functions of R. We denote the normalised time-like Killing vector field by — {e^'^, 0, 0, 0). Our 
solutions will be singularity and horizon free, so the coordinate system defined by eqn (||) is well behaved. We look 
for solutions of minimum energy. One can show [ p^f that the form of 'J compatible with this requirement is given by 



P 



exp 



an 



(9) 



where f2 is a real diniensionless constant and P is a real dimensionless function of R. A further requirement is that the 
eigenfunction P{R) possesses no nodes, so that the bosons are in their ground state. Then the product fla defines the 
ground-state energy of the bosons in the zero node eigenstate P{R)- From the line element (^) and the wave- function 
it is easy to show that T^u obeys the weak energy condition. 
Before writing the equations of motion, we define a new dimensionless radial coordinate 



aR 



(10) 



Then, using eqns (||) and 
system of coupled ODEs: 



1 



the independent components of the field and wave equations reduce to the following 

2 



20'r 



-ruj - 4(/)' - - (e^^ - l) 



(11) 



4(2cj + 3) 



g2Ap2 /^2g-^^ ^2uj + 5) + 2w - 1) P'^ {2uj + 1) 



+ 



4>''^ruj (jy'v'r 



,2A 



2r 



1) 



(12) 
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[e^^P^ (2 - n^e-^") + + + 0' fx' -y' - (13) 



2tj + 3 ^ ' V f 

and 
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P" = e'^P 1 - n'e-"") +P' [\' -y' --] , (14) 



r 



where the prime denotes 4-. To solve these equations we impose the boundary conditions 



e 



Ao _ 



1, P^ = 0, 0[, = O, Poo=0 (15) 



where throughout this paper the subscripts '0' and 'oo' denote values at r = and r — oo respectively. The first three 
conditions enforce regularity of the solutions at the origin, while the last condition ensures that the boson matter is 
localised and the solutions are asymptotically flat. For any given to the solutions may be parameterised by Pq and 
(j)oa- The field equations then become eigenvalue equations for fi. Note that the value of at infinity is arbitrary: 
eqns (|l^) to (^4|) are invariant under the rescaling v ^ v + k, SI —> Qe^ where k is constant. We use this freedom to 
set 

e"- = 1 (16) 

which fixes the scale of our time coordinate and consequently our unit of energy. The field equations automatically 
lead to the asymptotic conditions 

PL ^^'oo^ 4>'oo = Aoo - 0. (17) 

From eqns (H) and (^), is parallel to and has the non-zero component 

J' - F^e--. (18) 

Choosing a hypersurface orthogonal to and combining the above expression with eqn (|^) we have 

1 1"°° 

N^- e^-^nP'^r'^ dr (19) 

2 Jo 

where we have expressed A'' in units of M^i/a = fi/aG. Since there are no Maxwell terms present in the action (|^), 
the bosons have no electromagnetic charge and we interpret the quantity ^A'^ as the total boson particle number so 
that N is the total rest mass of the star. This latter quantity is simply the sum of the masses of the bosons that make 
up the star, as measured by some non-gravitational experiment. From the rest mass we define the Newtonian mass 

Mn ■■= Ne'l'°° (20) 

which measures the gravitational mass (or energy) of a star whose bosons are dispersed to infinity. 

Finally we note that to recover the GR equations of motion and their corresponding solutions we take the limit 
Lo oo which implies that (j) — > (jjcR, where (pcR is constant. Since we have included Newton's constant G in the 
action (|l|), we have used up the freedom to scale the dilaton (which is equivalent to rescaling the unit of mass) and 
so we have the additional requirement 4>gr — 0. 

III. QUASI-LOCAL AND ADM MASSES 

Operationally, in Newtonian theory one determines the active gravitational mass GM of a gravitating source by 
applying Kepler's third law. A test particle in a circular orbit of radius R and angular velocity dip/dt about the source 
measures an active mass 



GM = lim 

_R->oo 



R 



'dt 



(21) 



4 



This mass is a product of the gravitational couphng strength G and the rest mass M of the source. If the source is 
sphericaUy symmetric, eqn (|2^) is vahd for all R and it measures the gravitational mass GM{R) contained within the 
2-sphere of radius R. 

For metric theories of gravity the determination of mass is more difhcult. All we can discuss is some total gravita- 
tional mass, which is determined by the metric g^u and any additional tensor or scalar gravitational fields appearing in 
the theory. Furthermore, for metric theories other than GR, the Strong Equivalence Principle (SEP, defined in [|l^) is 
violated. Because of this, the mass we measure using Kepler's third law depends upon the test particle's gravitational 
binding energy. Hence, in an alternative theory of gravity we cannot in general separate the gravitational coupling 
strength from the rest mass of the source. 

Returning to the boson star solutions, we first discuss three definitions of quasi-local mass (QLM), adapted to the 
symmetries of the line element (^). As is well known, QLM definitions are of limited utility. One reason for this is 
that it is difficult to write a coordinate invariant conservation law for any quasi-local quantity. However, in a static 
spacetime, the time-like Killing vector field provides a natural splitting of the spacetime into a set of space-like 
hypersurfaces of constant time on which one can formulate an expression for the QLM. The existence of f guarantees 
that quantities defined on the constant time slices are conserved. There are still an infinite number of ways of defining 
QLM in a static spacetime. However, an acceptable definition of the mass M{U) of a sub-manifold U must satisfy the 
following three requirements: (a) M is non-negative, (b) M = everywhere in Minkowski spacetime, and (c) given 
two sub-manifolds U and V, such that U <ZV, then M{U) < M{V). We also feel that any global mass defined as the 
limit of some QLM is only valid if the QLM satisfies these conditions, and we use this criterion to select physically 
reasonable definitions of mass for the static boson star solutions we consider here. 

We first consider the generalised Schwarzschild mass, which is given by any one of the familiar definitions 



ms(r) = 




where ms{r) is measured in units of M^i/a and p is the total energy density defined by 

p G^,eC- (23) 

Equation ( p^ is the definition of mass adopted in all references on boson stars cited in the introduction except 
and it g ives an unambiguous definition of the energy of a spherically symmetric, asymptotically flat system in GR 
(see |16|] , p. 603): provided p is non- negative, nis satisfies conditions (a) to (c) outlined above and in the limit r — > oo, 
ms tends to the ADM mass. However, in Brans Dicke theory one or more of these conditions may be violated even 
for a physically reasonable energy momentum tensor. From eqns (^) and (^) the total energy density of the boson 
star is given by 

p = e-^^Goo = (l + I) e-'V" - + y {n^P^e-^^ + P^ + P'^e'^^) (24) 

where 

^-^^^-^^l^'^'V). (25) 

This expression contains non-positive definite terms so although T^j/ obeys the weak energy condition the total density 
may still be negative in some regions of a solution, leading to the condition m'g < (in violation of requirement (c)). 
This is indeed the case for the uj — —1 solutions discussed later, and there even exist weak field solutions for which 
ms < everywhere. The fact that the QLM ms may have a negative gradient leads one to conclude that the ADM 
mass, which as we shall see below is given by limfl^oo "m-s, istils to account for all of the energy of the spacetime. 

As an alternative to eqn (|2^), we may determine the gravitational mass enclosed within a two sphere S(r) by 
measuring the analogue of the gravitational acceleration of a test particle at radius r. More specifically, we consider 
circular geodesic orbits and use eqn ( |2l| ) to determine the mass enclosed by the orbit. This definition of mass is 
valid in a static, spherically symmetric spacetime. In this case we may project a particle orbit of constant r onto 
a hypersurface orthogonal to to form a circle about the centre of symmetry. The set of all orbits at this radius, 
projected in the same way, then form a closed two sphere S(r). The matter and gravitational fields contained within 
this sphere remain constant in time, and so the orbital mass is well defined. Since the SEP is violated in BD theory, 
the value of the orbital mass we measure depends upon the gravitational binding energy of the orbiting test particle. 

We consider two orbital masses: the Keplerian and the Tensor mass. A non-self gravitating test particle in a circular 
geodesic orbit in the geometry of eqn M) moves with an angular velocity 



5 



as measured by an observer at infinity. Use of eqn (^) tlien leads to the Keplerian mass function 

rriKir) = v'r^e^'' (27) 

in units of M'^J a. The Tensor mass we define similarly by considering the circular orbit of a test particle in the Einstein 
frame, whose metric g^i, is conformally related to the physical (Jordan frame) metric g^^, by the transformation 
g^iv = e~'^gfj,,y. The orbital angular velocity in the Einstein frame is given by 



'dt-'' \r{2-r^') ^^^> 
Kepler's third law (^ij) then gives the Tensor mass function 

2 — (p'r 

which may be interpreted as the mass measured by an orbiting test black hole in the Jordan frame As with the 
previous mass functions, the Tensor mass is given in units of M^j/a. Note that the conformal transformation used to 
define the Tensor mass is merely a computational device. We use it to find a frame in which the test particle's motion 
is geodesic and use this motion to define the Tensor mass. We could equally well define a different orbital mass by 
transforming to a different conformal frame. 

Numerical calculations for the u = — I BD solutions show that both niK and are non-decreasing functions of r. 
From eqns (^^ and (^8|), both mass functions are also non- negative definite (as a corollary eqn (|2^) implies that e"^ 
is a non-decreasing function of r). In Minkowski spacetime one has mx = mT = everywhere. Hence the Keplerian 
and Tensor masses satisfy conditions (a) and (b) given above and, at least for the a; = — 1 solutions studied later, 
they appear to satisfy condition (c). 

Lee [|8|] has shown that one can derive ADM-like masses in Brans-Dicke gravity from the superpotential 

^m[^"1/3 = (-g)e"('^~-'^) {g^"'g°''^ - g^'^'g"'^) (30) 

where n is usually taken to be an integer. Note that this expression is not unique due to the ambiguity in the choice 
of pseudo-tensor for the gravitational field. Furthermore, the index n is undetermined, which reflects the additional 
ambiguity we have in determining the contribution the dilaton field makes to the total gravitational energy momentum. 

The potential satisfies the conservation law H^'^"^ — and on integrating this expression and using Stoke's 
theorem twice we have the generalised ADM mass 

1 



M„ = 



where the integral is performed over the 2-sphere S(r) in the limit r oo and the metric must be expressed in 
rectangular (Minkowski) coordinates. Our definition of H^'^"^ differs from the one given in |^ by a factor of e'^'^°° 
which is included here so that M„ is expressed in units of M'^Ja. For the metric defined by eqn (||), the integral ( |3l| ) 
evaluates to 



ivin = um r/' I 1 — t 

r — >oo 

where 



M„= ^lim [.(l-e-^)]+^ (32) 



(1)1 = lim (r^0')- (33) 

r — >oo 

In the following Section we shall relate this expression to the asymptotic limits of the quasi-local masses. 
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IV. THE ASYMPTOTIC FORM OF THE SOLUTIONS 



Expanding the metric and mass functions in powers of 1/r about r — oo, one can show that to order 1/r the hne 
element is given by 

ds' = -(l- ^) dt^+(l + ^(MH^M) dr^ + r^dO^ + sin' 9 d^\ (34) 



where Mx is the hmit of eqn ( |27| ) as r — + cx) and 0i is defined in cqn (p3D. Since Um,,^o(^ v^') = this quantity may 
also be written as 

Equation ( ]34| ) is a special case of the vacuum BD solution given in and the asymptotic solution is determined 
by two parameters which we take here to be and Mk- For this form of the metric, the boson wave equation jl^ ) 
has the asymptotic solution 



P = r\- 

where 



l + O 



(36) 



2 



h = -\ - n{MK - 4>i) ^ — Mk, k= \/1-172. (37) 

The boson field falls off exponentially with r, far more rapidly than any other field. This justifies the use of the name 
"star" in describing these solutions: although the object has no well defined surface, the boson matter is still highly 
localised. The parameter k may be interpreted as the reciprocal radius of the star and must be real for an acceptable 
solution. This implies that < 1. 

These solutions are asymptotically flat so that for large r we have 1 — e'^^ ~ ^(1 — e^^^). Hence, setting n = in 
eqn (32) and comparing the result with cqn (p2|) gives 



Mo = Madai = lim ms (38) 

r — >QO 

and we see that, as mentioned above, the generalised Schwarzschild mass tends to the ADM mass in the Jordan frame. 
Taking the leading terms in the expansion of the Keplerian and Tensor masses, one can show that 

Mk = M4 = Madm + Mt = M2 = Madm + y . (39) 

Note that Mt is the ADM mass of the star in the Einstein frame, which is not equal to Madm since the derivative 
of the conformal factor e"*^ which relates the two frames is non- vanishing. Comer and Shinkai use the correct 
definition of mass (the Schwarzschild mass in the Einstein frame) but, as we have seen, this quantity is not the same 
as the Jordan frame ADM mass (this contradicts a statement given in |12|). Numerical studies show that 0' > for 
all w = — 1 solutions investigated. This implies that 

Madm < Mt < Mk (40) 

for these solutions. 

For each value of the mass M„ there is an associated fractional binding energy En per unit boson mass given by 

En - (Af„ - Mn)/Mn. (41) 

This expression is independent of the unit of mass chosen. The choice of which En correctly measures the binding 
energy of the star depends upon which mass value M„ corresponds to the true energy of the star. 



7 



V. DILATON RESCALING 



As mentioned above, the solutions may be parameterised by the pair Pq, (j)oo- The former quantity determines the 
energy density of the boson field at the origin and for each pair of parameters (Pq, 4'oo) there is a unique value of pair 



of He that gives a zero node solution with the boundary conditions (15^ and (16). The product Qe measures 



the bosons' ground state energy per unit boson mass at the origin, as seen by an observer at r = cx). This quantity 
increases monotonically with Pq and the pair (fJe""^" , 0oo) serves as an alternative parameterisation of the solutions. 
Note that fle~'^° is invariant under the rescaling of e'^ described before eqn (M). 

The equations of motion and the mass equations (^0|), (22), (|27|), (^) and (|32j) are invariant under the rescaling 



P 



P 
1 



where /c is a constant, while N re-scales under eqn 



as 



N 



N 



(42) 



(43) 



This rescaling swaps energy between the dilaton part of the gravitational field and the boson energy density is such 
a way as to leave the total gravitating energy invariant. Given a set of solutions 



S{Po;cl>oo)^{M{Po;cj)oo), W(Po;0oo)} 



(44) 



parameterised by Po for some fixed (poo, where M is any asymptotic or quasi-local mass value, eqn (^2|) generates a 
new physically distinct set 



5(Po;<^oo + 21ogfc) = <; Af(fcPo;<^oo), -^iV(fcPo; 0oo) 



(45) 



This new set of solutions have exactly the same mass and binding energy as the first set, but consists of stars of 
different total particle numbers with compensating changes in their gravitational coupling strength. Hence, eqn ( |45| ) 
may be used to compare boson star solutions in cosmological settings with different values of 0oo ■ 



VI. WEAK FIELD SOLUTIONS 



In this section we consider the behaviour of the zero node solutions in the weak field limit, which we define as the 
limit in which P{r) is small but non-zero. Following Kaup pO|, we write the boson field amplitude as 



P = en(a) 



(46) 



where < e ^ 1 and a is a new radial coordinate defined by a = ^/er. The function 11 is defined to have the 
boundary value Hq = 1, so that e measures the degree to which the solutions differ from flat spacetime. We use e as 
an expansion parameter and to first order in e we expand the metric functions, dilaton field and energy eigenvalue 
about flat spacetime as 



where r(Po,0oo) is positive and constant for each solution. Equations ( |ll| ) to (14) then reduce to 

2{uj + 2) 



2uj 



(47) 



(48) 



2ti> + 3 



(49) 



2uj + i 



0{e) 



(50) 
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and 



a^U')' = a^n{B + T) + 0{e) (51) 



where throughout this section a prime denotes 4- 



da ' 



Using eqns (19), (M), (|22|) and (M), the derivatives of eqns {U7h and (|29|), and the weak field equations of motion 



the quantities Mn, Madm, Mt and Mk may be written as the integrals 







2 



Jo 2uj + 3 



(53) 



Mt^M2^V^ da + Oie'^') (54) 

7o 2 

Afx = M4 = V^ r a'tfe*- (f^ii^ ^ 0(e^/^). (55) 
Jo 2w + 3 

Equation (^) shows that, in the weak field limit, the ADM mass is negative for lo < —1 — e. This is a consequence 
of treating </> solely as a matter field. The energy density may be written partly in terms of the weak field variables 
to give 

p = -^{a'<i>'y + 87rGe'f'T^,CC, (56) 

where we have shown the dilaton contribution to lowest non-zero order in e and one can show that the second term 
on the right is also of order e^. Equation ( pO| ) implies that (a^$')' is positive, so that the dilaton contribution to 
the total density in eqn (^) is negative. As uj decreases below w ~ — 1, the derivatives of the dilaton field increase 
in magnitude and the total density becomes dominated by the (negative) dilaton term. This leads to a negative 
Schwarzschild mass gradient for all r and so gives a negative quasi-local mass ms, and hence negative Madm- When 
UJ > —1 + e the coupling between the dilaton and the curvature is weak enough for the boson field density to dominate 
over the dilaton density, so that in this case both ms and Madm are positive. For — 1| ~ e we need to go to higher 
order in e to calculate Madm- However, numerical calculations show that AIadm increases smoothly from negative 
to positive values as lo increases over the interval (— 1 — e, — 1 -f e). 

From eqns (^2|), ( ^4|) and (|5^) we see that, asymptotically, the weak field solution describes the gravitational field 
about a central source of rest mass N so that to lowest order in e the Tensor and Keplerian masses may be written 
as the products M2 = G2N, M4 = G4N where the couphng strengths G2, G4 are given by 

G2=e^- G4 = e*~|^. (57) 

Hence the weak field solutions are Newtonian in the sense that the rest mass of the star may be separated from the 
gravitational coupling strength. However, this coupling is not unique and its strength depends upon the properties 
of the orbiting test particle that we use to measure the star's mass. The quantity G4, often quoted as the effective 
gravitational constant in the weak field limit (see for example |l^), is the coupling strength one would measure in 
a Cavendish experiment using non-self gravitating test particles. Performing the same experiment using test black 
holes, one would measure the coupling strength G2. 

We use eqns (^), (|5^) and (^5|) to calculate the fractional binding energies associated with the ADM and Keplerian 
masses. The result is 

Both of these quantities are independent of e in their leading term which is physically unreasonable: as £ ^ and the 
solution approaches flat spacetime, we would expect the fractional binding energy to vanish as it does in the GR case 
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pof . The fact that Eg and E4 behave in a non-physical way indicates that neither the ADM mass nor the Keplerian 
mass correctly describes the total energy of the star. 

As is apparent when considering eqns ( |5^ ) and (jsj), to calculate E2 we need to express Mn and Mt to order e^^^. 
It turns out that the ©(e^^^) terms of both of these quantities depend only upon the 0{e) weak field variables defined 
m eqn ®. One can then show that the fractional binding energy associated with the Tensor mass is given by 



E2=-—r 



M Jo 2w + 3 



$'(w + 2) + B' 



,(2t^ + l) 



da + 0(6"^), 



where we have defined 



da 



(59) 



(60) 



so that y/eAf is the rest mass to lowest order in e. We need to simplify eqn i^^. Integrating eqns (^) and ( |50| ) one 
can show that 



B' 



2{uj + 2) 



and substituting this result into eqn (p9|) gives 



El — — —r 



°° a3n2 



B' da + 0{e^). 



(61) 



(62) 



Integrating both a^H'^ and all' (a^II')' by parts and using eqn ( ^l|) one can show that 

poo pOO 

/ a^U^B' da = -2 / a^n^{T + B) da. 
Jo Jo 

by parts and using eqn (^8| ) one can show that 



(63) 



Similarly, integrating both a B and 



aB' ia^B') 



a^U^B da = -2 / a^H^B' da 



Substituting eqn (64) into (|63[) gives 



2r 



2r 



a^n^B' da= — a^U^ da = —N 



3 Jo 

where we have used eqn (|o|) to obtain the last equality. Combining the above result with eqn ( |6^ ) gives 

E,^-'-^+0{e^). 



(64) 



(65) 



(66) 



This is identical to the expression for the GR binding energy found by Kaup |pO| and, in particular, is independent 
of u) and vanishes as e ^ 0. Hence all weak field BD boson stars are energetically stable. 



VII. STRONG FIELD SOLUTIONS 



We have numerically integrated the field equations (|llj) to (14) for the case w = — 1 over the parameter range 
< Po < 1.6. Results of the integration are given in Figure 1, which shows the behaviour of the three asymptotic 
masses we have discussed and how they relate to the Newtonian mass Mm- We have chosen the boundary value 
g<^oo = 1 go that Mm is numerically equal to N and our results may may be compared with those of other authors, 
although this choice is not physically realistic: by setting lo = —1 we are implicitly describing boson star solutions 
in the early (string-theory) Universe. In general, ST cosmological solutions show that the cosmological value of 
decreases with cosmic time, so that we should choose the boundary value 6*^°° 3> 1. However, from the scaling relation 
(O) , changing the asymptotic value of cf) merely changes the scale of the the horizontal axis of the figure (although 
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for 6*^°° ^ 1, Mat and iV are no longer equal). Note also that the fractional binding energy of any solution is invariant 
under the rescaling of (^^o (this result is true for all values of w). 

The four mass curves exhibit qualitatively similar behaviour to the ADM mass and particle number curves found 
for GR boson stars. Each varies smoothly with Pq and each mass curve reaches an initial maximum before falling off 
with increasing Pq and tending towards low amplitude oscillations about some low mass value. There are, however, 
several quantitative differences between the GR and BD solutions. 

All of the solutions investigated possess a region in which —0" is large and negative and dominates over the boson 
field contribution to the total energy density. This implies that p < in this region, which leads to the condition 



m'g < there, in violation of condition (c) given in Section III. The quasi- local Tensor and Keplerian masses niT and 



•niK were found to be increasing functions of r for all solutions. 

An obvious feature of the Figure is that the asymptotic masses M^/^m , and Mj- differ remarkably from each 
other. This is a result of choosing such a low value of w, which means that the gradient of is large, even in the 
asymptotic limit. Gonsequently, for these solutions 0i is of the same order of magnitude as Madm and from eqn ( |39| ) 
the difference between the values of the three masses at each Pq is large. As lo increases, the dilaton field becomes 
less strongly coupled to the curvature and the ratio (Pi/Madm decreases for a given N, which implies that both Mk 



and Mt tend towards Madm- From eqns ( |53| ) and (54), for a given Pq the fractional difference between Madm and 
Mt in the weak field limit is given by 

Mt - Madm ^ 1 

Mt ^ 2uj + 3' ^ ^ 

Further numerical work suggests that this relation is approximately true for strong field solutions, with an accuracy 
that increases with oj. For the oj — 6 coupling chosen by Gunderson & Jensen H], the fractional difference between 
Madm and Mt is around 1/15, which implies that the boson star masses that they quote are about 93% of their true 
values. For the choice w = 500 implied by the current observational constraints, the fractional difference between the 
two masses is less that 10~^ which, as pointed out in ||l^, is negligible. 

A second obvious feature of the Figure is that the Madm, Mk and Mt curves have extrema located at different 
values of Pq . This feature persists for all finite uj and is only absent in the limit w — > cxd where the three masses are 
identical. However, it does appear from the figure that the extrema of Mt and Mn do occur at the same values of 
Pq, which suggests that Mt is the correct definition of the total energy of the star. In the following Section we give 
a proof that this is indeed the case for all values of oj. 



VIII. COINCIDENCE OF EXTREMAL TENSOR MASS AND EXTREMAL PARTICLE NUMBER 

SOLUTIONS 



We consider variations SMt and SMn in the Tensor and Newtonian masses induced by the parameter change 
Pq ^ Pq + SPq with (poo held fixed. Since each pair Pq, 4>oo gives a unique asymptotically flat zero node solution to 
the equations of structure (|ll|) to (Q), we may equivalently consider arbitrary variations in Mt and M^ subject to 
the constraint that the field equations and the boundary conditions (15) to ( p7| ) hold. We shall then end up with a 
simple relationship between SMt and SMn- Note that because of our choice of coordinates, the variation S commutes 
with both integration and differentiation with respect to r. 

We start by taking the variation in the rest mass N. From the definition (n3) we have 



6N = 



Using a similar method to the one outlined in 
way. One can easily show that 



/ —[e^-''np^6X + e^6{e-''nP^)] dr. (68) 
Jo 2 

we rewrite the second term in the above integrand in the following 



np2 



-"npsp. 



(69) 



We solve eqn (24) for the term proportional to il^P^ and substitute the result into the right hand side of the above 
equation. Then, substituting this result into eqn (p8[) gives 



6N 



r 

2n 



e-'+^-'f' [5p + e^-2A (g2A-2.^2p2 ^ p,2 + (2 + ^) e-^0'2) S\ + e^p (^e'^-fl^ - l) 6P 



-2ur,2 n2 



n^p^ + P^ 



-2\pl2 



~2\ r,l 



P'SP' - (2 + uj)e-^^(l)'S(l)' + (5(A(/)) 



dr. 



(70) 
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Integrating the 5P' term by parts and using the wave equation (|TJ) one obtains an expression that cancels the 5P 
term in eqn ([70|). Integrating the <5(A0) term by parts we have 



2n 



Q Zi ' 



- lim 

r — >oo 



2n 



Substituting this result into eqn ( |70| ) gives 

rOG 2 

SN^ ^e''+^-^ [Sp + (e'^^il^P^ + e-^^P'^ + (1 + Lu)e-^^-^^'^ - e'^Acf) + e'^^-'^i^' (f>') S\ 
Jo 2i2 



+ lim 

r~*QO 



2n 



Integrating the Scj)' term by parts and substituting the result back into eqn (|7|) we have 
Jo 2f2 



2n 



-S (e-V20') 



(71) 



(72) 



(73) 



where we have used eqn ( p4[ ) to rewrite part of the coefficient of S(j) in terms of the energy density p and rewritten 
some of the derivatives as the wave operator □. This latter quantity has the explicit form 



for any function f{r). 

From the definition of the Einstein tensor we have 



(74) 



(75) 



where Rap is the Ricci tensor. Rearranging this equation, substituting in the explicit form of Rai3^°'£,'^ = e ^'^Rqq 
and using the definition (E3h wc have 



n 

2" 



(76) 



where we have used eqn ( |7^ ) to obtain the second equality. Solving eqn (Q) for TZ, substituting the result into the 
above equation and evaluating the a4> term gives 



p + ^uoer'^^cf)''^ -Uv - ujUcf) = 0. 

Hence the Scj) term in eq n ([73| ) vanishes. 

Combining eqns (O), (p2p7(p3h and (p4) one can show that 



„2A-2.^2p2 ^ p,2 g2A-0^^ ^ ^-0^, y ^ ^ ^ — (5-+^- 

r dr ^ 



Substituting this result into the remaining part of eqn (73) gives 

SX 



SN 



2n 



e^'-^-'t'Sp 



.e-2A:^/g.+A- 



dr 



dr + lim 

r — *oo 



2rj 



Taking the variation of eqns ( p^ ) one can show that 



(77) 



(78) 



(79) 



(80) 
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and 



„2A 



i5A = 5mc. 



Substituting both of these results into eqn ([79^ gives 



5N = 



1 d 



ndr 



5ms) dr + lim 



2Vl 



Evaluating the integral and using the boundary conditions (h5|) to (nju we have 



SN : 



lim 

r — >oo 



5ms 5{r'^4>')'^ 



1 



(81) 



(82) 



(83) 



where we have used eqns (|3^), ( |38D and (|39|) to obtain the second equality. Finally, using the definition (^0|) we have 

SMt = n SMn- (84) 

The above equation gives the relationship between the variations in the Tensor and Newtonian masses under the 
parameter change Pq Pq + SPq with held fixed. Taking the limit SPq — > we have 



dMr 



dPn 



dJVlN 



dPo 



(85) 



Hence for any uj and (j)oo, a solution that extremises the Newtonian mass is also one of extremal Tensor mass. In the 
GR limit, Mt Madm and we recover the result derived in B]. 



IX. CONCLUSIONS 



We have examined three alternative definitions of mass in ST gravity: the ADM mass Madm, the Keplerian mass 
Mk and the Tensor mass Mt- In a static, spherically symmetric spacetime, all three may be written as the asymptotic 
limit of some quasi-local mass defined on the hypersurfaces of constant time. In the case of the ADM mass we have 
written this as the asymptotic limit of the Schwarzschild mass ms which is defined in the usual way as the integral 
of the total energy density G^i^^'*^'^ over a bounded subset of a constant time hypersurface. For a static spacetime 
in GR, the Schwarzschild mass defined in this way gives a perfectly satisfactory notion of quasi-local mass. We have 
analysed the behaviour of the asymptotic and quasi- local masses for BD boson stars in the strong coupling (co — — 1) 
case and found that both ms and Madm have undesirable physical properties even though the boson field energy 
momentum tensor obeys the weak energy condition. In addition, for strong couplings, Madm < Mn for all values 
of the central boson field amplitude. If one assumes that the energy of the star is given by Madm, one would have 
to conclude that all stars in the strong coupling case are energetically stable, regardless of their central density. This 
contrasts strongly with the behaviour of boson star solutions in both GR and in weakly coupled BD theory. In all of 
these cases, one finds that the solutions become unstable above a certain value of the central density. These results 
imply that Madm does not give a correct description of the energy of the spacetime, primarily because it does not 
correctly include the dilaton field's contribution to the total energy. The dilaton, being an extra component of the 
gravitational field, must be treated on a slightly different footing to the normal matter fields. 

Both Mt and Mk adequately include the contribution made by the dilaton to the active gravitational mass of the 
boson star as seen by two different kinds of orbiting test particles. However, as pointed out by Lee ||], out of the 
range of possible masses one can derive from the superpotential (pO|), it is only the Tensor mass that is conserved in a 
non-static isolated source. This implies that the Tensor mass may be the true physical energy of the star. However, 
the analysis given by Lee does not show that Mt uniquely satisfies a conservation law: one can conceive of other 
definitions of energy that are also conserved. In this work we have shown that, for a one parameter set of boson star 
solutions in BD theory, Mt is unique in that it is the only conserved quantity whose extremal points coincide with 
those of the boson particle number. In this respect, Mt is the analogue of the ADM energy used in GR and this 
behaviour supports the conjecture that Mt uniquely describes the energy of the system. We see no reason why these 
results should not hold for stellar objects composed of matter other than bosons, or in more general scalar-tensor 
theories. 

For the weak coupling case, the numerical difference between Mt and Madm is small, so in practice there is little 
accuracy lost in identifying the energy of a boson star with Madm- However, in the early Universe and in more 
general scalar-tensor theories, one may have to consider strong couplings and in these cases it is important that the 
correct definition of mass is used. 
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FIG. 1. Mass curves for the oj = —1 solutions with e'*°° = 1. The curves are parameterised by Po, the amplitude of the 
boson wave function at the origin. Four mass curves are shown: the Newtonian mass Mn (solid line), the ADM mass Madm 
(dashed line), the Tensor mass Mr (dotted line) and the Keplerian mass Mk (dot-dashed line). All masses are expressed in 
units of Mpi/a. 
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